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The self-consistent theory of Anderson localization of quantum particles or classi- 
cal waves in disordered media is reviewed. After presenting the basic concepts of 
the theory of Anderson localization in the case of electrons in disordered solids, 
the regimes of weak and strong localization are discussed. Then the scaling theory 
of the Anderson localization transition is reviewed. The renormalization group 
theory is introduced and results and consequences are presented. It is shown 
how scale-dependent terms in the renormalized perturbation theory of the in- 
verse diffusion coefficient lead in a natural way to a self-consistent equation for 
the diffusion coefficient. The latter accounts quantitatively for the static and 
dynamic transport properties except for a region near the critical point. Several 
recent applications and extensions of the self-consistent theory, in particular for 
classical waves, are discussed. 



1.1. Introduction to Anderson localization 
1.1.1. Brief historical review 

The localization of quantum particles by a static random potential, or of classical 
waves by random fluctuations of the medium, is one of the most intriguing phe- 
nomena in statistical physics. The key ingredient of localization, wave interference, 
was introduced in P. W. Anderson's seminal paper "Absence of diffusion in certain 
random lattices" } There it was shown that electrons may be localized by a random 
potential, so that diffusion is suppressed, even in a situation where classical particles 
would be dclocalizcd. The fundamental reason for the localizing effect of a random 
potential on quantum particles or classical waves is the multiple interference of wave 
components scattered by randomly positioned scattering centers. The interference 
effect takes place, as long as the propagation is coherent. 

It is interesting to note that the first application of the idea of localization 
concerned the spin diffusion D of electrons and not the electrical conductivity a. 
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Anderson considered a tight-binding model of electrons on a crystal lattice, with 
energy levels at each site chosen from a random distribution. 1 The traditional view 
had been, that scattering by the random potential causes the Bloch waves to lose 
well-defined momentum on the length scale of the mean-free path £. Nevertheless, 
the wavefunction was thought to remain extended throughout the sample. Anderson 
pointed out that if the disorder is sufficiently strong, the particles may become 
localized, in that the envelope of the wave function ip(r) decays exponentially from 
some point Vq in space: 

| V(r) |~ exp(| r - r | /£), (LI) 

where £ is the localization length. 

There exist a number of review articles on the Anderson localization problem. 
The most complete account of the early work was presented by Lee and Ramakr- 
ishnan. 2 The seminal early work on interaction affects is presented in Ref. 3 A 
complete account of the early numerical work can be found in Ref. 4 A path integral 
formulation of weak localization is presented in Ref. 5 Several more review articles 
and books are cited along the way. In the following we will use units with Planck's 
constant h and Boltzmann's constant fcs equal to unity, unless stated otherwise. 

1.1.2. Electrons and classical waves in disordered systems 

The wavefunction tp(r) of a single electron of mass m in a random potential V(r) 
obeys the stationary Schrodingcr equation 

(-|^V 2 + V(r)-i?)^(r) = 0. (1.2) 

In the simplest case V(r) may be assumed to obey Gaussian statistics with 
(V(r)V(r')) = (V 2 )5(r — r') 1 but many of the results presented below are valid 
for a much wider class of models. Electrons propagating in the random potential 
V(r) will be scattered on average after a time r. For weak random potential the 
scattering rate is given by 

- = nN(E) (V 2 ) (1.3) 

T 

where N(E) is the density of states at the energy E of the electron. In a metal 
the electrons carrying the charge current are those at the Fermi energy E = Ep. 
Within the time r the electron travels a distance I = VpT, where vp is its velocity. 

In close analogy the wave amplitude ip( r ) of a classical monochromatic wave of 
frequency ui obeys the wave equation 

(c^) +V2 V (r)= °- (L4) 

Here c(r) is the wave velocity at position r in an inhomogeneous medium, assumed 
to be a randomly fluctuating quantity. The main difference between the Schrodingcr 
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equation and the wave equation is that in the wave equation the "random potential" 
l/c 2 (r) is multiplied by uj 2 , so that disorder is suppressed in the limit ui — > 0. By 
contrast, in the quantum case disorder will be dominant in the limit of low energy 
E. A further difference may arise if the wave amplitude is a vector quantity as, e.g., 
in the case of electromagnetic waves in d — 3 dimensions. 

In real systems particles or wave packets are not independent, but interact. 
Electrons are coupled by the Coulomb interaction, leading to important effects that 
go much beyond the single particle model. Similarly, wave packets interact via 
nonlinear polarization of the medium. Apart from these complications, the physics 
of electronic wave packets and classical wave packets is quite similar. In the following 
we will present most of the discussion in the language of electronic wave packets. 

1.1.3. Weak localization 

The all-important effect of wave interference is most clearly seen in the limit of weak 
scattering, where it already may cause localization, but only in reduced dimensions. 
While it is difficult to observe full localization at finite temperature T, on account 
of the effect of interactions limiting the phase coherence, the dramatic signatures 
of localization are visible at finite T in the form of "weak localization" , 6,7 

An electron or a wave packet moving through a disordered medium will be 
scattered by the random potential on the average after propagating a distance £, 
the mean-free path. On larger length scales the propagation is diffusive. Weak 
localization is a consequence of destructive interference of two wave components 
starting at some point and returning to the same point after traversing time-reversed 
paths. Let the probability amplitudes for the wave packet to move from point ro 
along some path C\ back to ro be A\ and along a different path C2 be A2, then the 
transition probability for the particle to move either along C\ or along C2 will be 

w =| At + A 2 | 2 = w ci + w int , (1.5) 

where w c i =| A\ | 2 + | Ai | 2 and wt n t = 2Re(A\A2) . For any two paths the inter- 
ference term Wi n t may be positive or negative, and thus averages to zero. However, 
if A2 = A r is the amplitude of the time-reverse of path A\ = A and if time reversal 
holds, then A = A r , i.e., the probability of return w is enhanced by a factor of two 
compared to the probability w c i of a classical system: 

w = 4 I A | 2 = 2w d . (1.6) 

In that case the probability for transmission is reduced, which leads to a reduced 
diffusion coefficient and a reduced conductivity. One may estimate the correction 
to the conductivity in the following qualitative way. The relative change of the 
conductivity a by the above interference effect is equal to the probability of inter- 
ference of two wave packets of extension A, the wavelength, after returning to the 
starting point. The infinitesimal probability of return to the origin in time t of a 
particle diffusing in d dimension is given by {AnDt)~ d / 2 d z r where D is the diffusion 
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coefficient. Since the volume of interference in the time interval [t, t + dt] is X d ~ 1 vdt, 
where v is the velocity of the wave packet, one finds the quantum correction to the 
conductivity 8 a as 6 ' 7 



6a r* v\ d ~ l dt 



-c 2 j In^/r), d = 2 (1.7) 
-ci(J$-l), d=l. 



a 7 T {AnDt) d / 2 

Here -D = ^i> 2 r is the diffusion constant, £ — vt is the mean-free path, r is the mean 
time between successive elastic collisions, ao = e 2 nT /m is the Drude conductivity 
with n as the particle density, and Ci are constants of order unity. The upper limit 
of the integral is the phase relaxation time r^, i.e., the average time after which 
phase coherence is lost due to inelastic or other phase-shifting processes. For weak 
localization processes to exist at all, the inequality ^> t must hold. We note that 
the correction in three and two dimensions depends on the ratio of wavelength A to 
mean-free path £, and gets smaller in the limit of weak disorder, where \jl -C 1. In 
two and one dimension the correction grows large in the limit t/t^ — > since one 
expects the phase relaxation rate 1 /t$ for a system in thermodynamic equilibrium 
to go to zero for T — > 0. By contrast, in some cases a plateau behavior of l/r^ 
as a function of temperature has been found experimentally, which gave rise to the 
speculation that the zero point fluctuations may cause decoherence. However, given 
a unique ground state, it is difficult to understand how a particle in the system may 
loose its phase coherence. Several physical mechanisms that may lead to a plateau 
of 1/t<£ have been identified. For a recent discussion of these issues see Rcf. 8 

With t/t^ — >• for T — >• the weak localization quantum correction will be 
large in any system in d = 1,2, no matter how weak the disorder. As we will 
see, this behavior signals the fact that there are no extended states in d = 1,2 
dimensions. The characteristic length L$ over which a wave packet retains phase 
coherence is related to by the diffusion coefficient = -J 'Dt$. In systems of 
restricted dimension, e.g., films of thickness a or wires of diameter a, the effective 
dimensionality of the system with respect to localization is determined by the ratio 
Lfft/a. Namely, for L$ <C a the system is three-dimensional (3c?), while for > a 
diffusion over time takes place in the restricted geometry of the film or wire, and 
the effective dimension is therefore 2 or 1. 

1.1.4. Strong localization and the Anderson transition 

The appearance of localized states is easily understood in the limit of very strong 
disorder: localized orbitals will then exist at positions where the random potential 
forms a deep well. The admixture of adjacent orbitals by the hopping amplitudes 
will only cause a perturbation that does not delocalize the particle. The reason 
for this is that nearby orbitals will have sufficiently different energies so that the 
amount of admixture is small. On the other hand, orbitals close in energy will in 
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general be spatially far apart, so that their overlap is exponentially small. Thus, we 
can expect the wave functions in strongly disordered systems to be exponentially 
localized. Whether the particles become delocalized when the disorder strength 
is reduced, is a much more complex question. In one dimension it can be shown 
rigorously that all states are localized, no matter how weak the disorder. 9-11 

In three dimensions, the accepted view is that the particles are delocalized for 
weak disorder. In general, localized and extended states of the same energy do 
not coexist, since in a typical situation any small perturbation would lead to hy- 
bridization and thus to the derealization of a localized state. We can therefore 
assume that the localized and extended states of a given energy are separated. For 
increasing disorder strength rj there will then be a sharp transition from delocalized 
to localized states at a critical disorder strength rj c . A qualitative criterion as to 
when an Anderson transition is expected in 3d-systems has been proposed by Ioffc 
and Regel. 12 It states that as the mean free path I becomes shorter with increasing 
disorder, the Anderson transition occurs when I is of the order of the wavelength 
A of the particle (which amounts to the condition hpl ~ 1 in metals, where hp is 
the Fermi wave number). As we will see later, in Id or 2d systems £ may be much 
longer than the wavelength and the particles are nonetheless localized. In fact, the 
relevant mean free path here is the one with respect to momentum transfer. A 
similar situation exists when we fix the disorder strength, but vary the energy E. 
Electrons in states near the bottom of the energy band are expected to be localized 
even by a weakly disordered potential, whereas electrons in states near the band 
center (in d = 3) will be delocalized, provided the disorder is not too strong. Thus 
there exists a critical energy E c separating localized from delocalized states, the 
so-called mobility edge. 13 ' 14 The electron mobility as a function of energy is iden- 
tically zero on the localized side (at zero temperature), and increases continuously 
with energy separation | E — E c | in the delocalized, or metallic, phase. The con- 
tinuous character of this quantum phase transition, termed Anderson transition, is 
a consequence of the scaling theory to be presented below. 

Historically the continuous nature of the metal-insulator transition in disordered 
solids has been a point of controversy for many years. According to an earlier theory 
by Mott 13,14 the conductivity changes discontinuously at the transition, such that 
a "minimum metallic conductivity" exists on the metallic side of the transition. 
Numerical simulations have shown beyond doubt that the transition is instead 
continuous, at least in the absence of interactions. 

In the much more complex situation of interacting electrons one finds for 
the Hubbard model without disorder, using the Dynamical Mean-Field Theory 
(DMFT), that the Mott- Hubbard metal-insulator transition is discontinuous at fi- 
nite temperatures, and that it becomes continuous in the limit T — > 0. 15,16 For 
the Hubbard model in the presence of disorder ("Anderson-Hubbard model") at 
T = the situation is similar: the Mott-Hubbard metal-insulator transition is 
discontinuous for finite disorder and becomes continuous in the limit of vanishing 
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disorder. 17,18 

1.2. Fundamental theoretical concepts of Anderson localization 

The Anderson localization transition is a quantum phase transition, i.e., it is a tran- 
sition at zero temperature tuned by a control parameter, e.g., the disorder strength, 
particle energy, or wave frequency. Unlike other quantum phase transitions, the 
Anderson transition does not have an obvious order parameter. Nonetheless, there 
exists a dynamically generated length scale, the localization or correlation length £, 
which tends to infinity as the transition is approached. Therefore, by drawing an 
analogy with magnetic phase transitions, Wcgner early on proposed scaling proper- 
ties. 19 Later, he formulated a field-theoretic description of the Anderson transition 
in the form of a non-linear sigma model (NLerM) of interacting matrices (rather 
than vectors, as for magnetic systems). 20 The NLcrM was later formulated in the 
mathematically more tractable supersymmetric form. 21 

1.2.1. Scaling theory of the conductance 

Wegner 19 argued that the Anderson localization transition should be described in 
the language of critical phenomena of continuous (quantum) phase transitions. This 
requires the assumption of a correlation length £ diverging as a function of disorder 
strength 77 at the critical point 

m h v - vc r ■ (i-8) 

The conductivity is then expected to obey the scaling law 

a(n)~e- d ~(Vc-ri) s ; V < Vc d > 2. (1.9) 

This follows from the fact that er, written in units of e 2 /(27rS), has dimension 
(l/length)^ 2 , and the only characteristic length near the transition is the correla- 
tion length £. By comparing the conductivity exponent s with the exponent of £ 
one finds 

s = v(d-2). (1.10) 

On the other hand, the conductance g of a d-dimensional cube of length L, which 
for a good metal of conductivity a is given by g(L) = aL d ~ 2 , must obey the scaling 
property 

g( V ;L) = <t>(L/0- (1.11) 



This means that g is a function of a single parameter L/£, so that each value of L/£ 
corresponds to a value g. 
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1.2.2. Renormalization group equation 

It then follows that g(L) obeys the renormalization group (RG) equation 

!&-«•>■ <" 2 > 

where j3(g) is a function of g only, and does not depend on disorder. In a landmark 
paper, Abrahams, Anderson, Licciardcllo and Ramakrishnan 6 proposed the above 
equation and calculated the /3-function in the limits of weak and strong disorder. 
A confirmation of the assumption of scaling was obtained from a calculation of the 
next-order term. 7 

At strong disorder one expects all states to be localized, with average localization 
length £. It then follows that g(L) is an exponentially decreasing function of L: 

g(L) ~ exp(-L/0. (1.13) 

In comparison with the ohmic dependence g ~ L d ~~ 2 this is a very non-ohmic be- 
havior. The ^-function is then given by 

(3(g) ~ \n(g/g c ) < 0. (1.14) 

At weak disorder one finds from g ~ L d ~ 2 that 

P(g) = d-2. (1.15) 

The important question of whether the system is delocalized (metal) or localized 
(insulator) may be answered by integrating the RG equation from some starting 
point Lq, where g{Lo) is known. Depending on whether /3(g) is positive or negative 
along the integration path, the conductance will scale to infinity or to zero, as L 
goes to infinity. 

In d = 3 dimensions one has (3(g) > at large g, but (3(g) < at small g. Thus, 
there exists a critical point at g = g c , where (3(g c ) = 0, separating localized and 
delocalized behavior. 

On the other hand, in d = 1 dimension one has (3(g) < at large and small g, 
and by interpolation also for intermediate values of g, so that there is no transition 
in this case and all states are localized. 

The dimension d = 2 apparently plays a special role, as in this case (3(g) — > 
for g — > oo. In order to determine whether (3 > or < for large g one has to 
calculate the scale dependent (i.e., i-dependent) corrections to the Drude result at 
large g. This is precisely the weak localization correction already mentioned above. 
For a system of finite length L < we should replace ^- = DL^ 2 in eq. (1.7) by 
DL~ 2 , leading to 

g(L)=a -a ln(^), (1.16) 

where a diagrammatic calculation 6 gives a = 2/tt and ctq = £/Xf (in units of e /H; 
Xf is the Fermi wave length) has been used. It follows that 

P(g) = -- g , d = 2, (1.17) 
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so that we can expect /3(g) < for all <?, implying that again all states are localized. 
This result is valid for the "usual" type of disorder, i.e., in case all symmetries, in 
particular time reversal symmetry (required for the weak localization correction to 
be present) are preserved. If time- reversal invariance is broken, e.g., by spin- flip 
scattering at magnetic impurities, the weak localization effect is somewhat reduced 
in dimensions d = 2 + e, e << 1, but is not completely removed. The first correction 
term in the j3- function is then proportional to — 1/g 2 (see, e.g., Ref. 21 ) implying 
that all states are still localized (in d = 3 dimensions the leading correction term 
is again ~ 1/g; see Ref. 22 ). In the presence of a magnetic field the situation is 
more complex, since the scaling of the Hall conductance is coupled to the scaling 
of g. As a result, one finds exactly one extended state per Landau energy level, 
which then gives rise to the quantum Hall effect. 23 On the other hand, if spin- 
rotation invariance is broken, but time-reversal invariance is preserved, as is the 
case of spin-orbit scattering, the correction term is proportional to +l/g, i.e., it is 
anti-localizing. In this case the /3-function in d = 2 dimensions has a zero, implying 
the existence of an Anderson transition. 24 

1.2.3. Critical exponents 

In the neighborhood of the critical point at g = g c in d = 3 we may expand the 
/^-function as 



f 



y 



P{3) = - 9 — jL > I 9 -9c |« 9c (1.18) 



9c 



Integrating the RG equation for g > g c from g(l) = go to — > I at large L we find 
g(L) = crL, where 

<r~j(g(l)-g c y. (1.19) 

Since [g(€) — g c ] cx (r) c — J]), we conclude that the inverse of the slope of the /?- 
function, y, is equal to the conductivity exponent s = y. 
Similarly, one finds on the localized side (g < g c ) 

g(L) ~ .g c exp [ - c(g c - g(£)) v L/i\ ~ g e exp(-L/S), (1.20) 

from which the localization length follows as 

Z~e\v-vc\- y . (i.2i) 

The critical exponent v governing the localization length is therefore v = y = s in 
d = 3 dimensions. 

Since the critical conductance g c = 0(1) in d = 3, there exist no analytical 
methods to calculate the /3-function in the critical region in a quantitative way. 
A perturbative expansion in 2 + e dimensions, where g c 3> 1, is possible, but the 
expansion in e is not well-behaved, so that it cannot be used to obtain quantitative 
results for s and v in d = 3. There exist, however, reliable results on v from 
numerical studies, according to which s = v = 1.58 ± 0.02. 4 ' 25 
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1.2.4. Dynamical scaling 

The dynamical conductivity cr(cj), i.e., the a.c. conductivity at frequency uj, in the 
thermodynamic limit in d = 3 obeys the scaling law 26,2 

a(uj;r)) = ^(L u /0, (1.22) 

where the scaling function $ has been introduced in eq. (1.11). Here is the typical 
length which a wave packet travels in the time of one cycle, 1/uj. Since the motion 
is diffusive it obeys L w = ^D(uj)/uj. It is important to note that the diffusion 
coefficient D{uj) is energy scale dependent and is related to the conductivity via the 
Einstein relation 

cr(w) = hN{E)D{uj), (1.23) 

where N(E) is the density of states at the particle energy E. 

At the Anderson transition, where £ — > oo, we expect cr(uj) to be finite. It follows 
that lim^oo 3>(Z/ w /£) ~ £/ ' L u and consequently 

a(ui;n) ~ -!- , r) = Vc- (1-24) 
This is a self-consistent equation for u(uj), with solution 

ctH-w 1 / 3 , 77 = 770. (1-25) 

More precisely, in the above expressions uj should be replaced by the imaginary 
frequency —iiu, such that cr(w) is a complex- valued quantity. 

In a more general notation, introducing the dynamical critical exponent z by 
<j{lo) ~ w 1 / 2 , we conclude that z = 3. The dynamical scaling is valid in a wide 
neighborhood of the critical point, defined by uj > ^(£/£)~ 2 ~| 77 — r\ c \ vz , where 
i/z sa 4.8. This scaling regime is accessible in experiment, not only by measuring the 
dynamical conductivity directly, but also by observing that at finite temperature 
the scaling in uj is cut off by the phase relaxation rate l/r^. 27 Therefore, assuming 
a single temperature power law 1/r^ ~ T p , one finds the following scaling law for 
the temperature dependent d.c. conductivity 

cr(T; 77) - T p / 3 ^> T (£_T p / 3 ). (1.26) 

Using this scaling law one may in principle determine the critical exponent v from 
the temperature dependence of the conductivity in the vicinity of the critical point. 
In the case of disordered metals or semiconductors, where studies of this type have 
been performed, the effect of electron-electron interaction has to be taken into ac- 
count. One major modification in the above is that the Einstein relation is changed. 
Namely, the single-particle density of states (which is not critical) is replaced by 
the compressibility dn/d/j,, with 71 as the density and /x as the chemical potential, 
which in the presence of the long-range Coulomb interaction is expected to vanish 
at the transition, i.e., the system becomes incompressible. Another change is that 
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the frequency cutoff is given by the temperature. The critical exponents determined 
from experiment vary widely, from s = 0.5 (Ref. 28 ) and s = 1 (Ref. 29 ) to s = 1.6 
(Ref. 30,31 ), and from z = 2 (Ref. 30 ) to z = 2.94 (Ref. 31 ). 

1.3. Renormalized perturbation theory of quantum transport in dis- 
ordered media 

The field-theoretic description in terms of the nonlinear a model (NLcrM) mentioned 
in the beginning of Sec. 1.2 is believed to be an exact framework within which 
the critical properties of the Anderson transition may be, in principle, calculated 
exactly. The mapping of the initial microscopic model onto the NLcrM requires a 
number of simplifications, so that the noncritical properties like the critical disorder 
r) c , the behavior in anisotropic systems, or systems of finite extension arc no longer 
well represented by this model. In addition, it is not known how to solve the NLcrM 
in cases of major interest, such as in d = 3 dimensions. 

It is therefore useful to consider approximation schemes, which on one hand 
keep the information about the specific properties of the system and on the other 
hand account approximately for the critical properties at the transition. Such a 
scheme, the self-consistent theory of Anderson localization, is available at least for 
the orthogonal ensemble (in which both, time reversal and spin rotation symmetry 
are conserved). This approach has been developed by us in Rcfs. 32,33 and was 
reviewed in Ref. 34 It may be termed "self-consistent one-loop approximation" in 
the language of rcnormalization group theory but has, in fact, been derived following 
a somewhat different logic as will be discussed below. 

The appropriate language to formulate a microscopic theory of quantum trans- 
port or wave transport in disordered media is a renormalized perturbation theory in 
the disorder potential. The building blocks of this theory for the model defined by 
cq. (2) arc (i) the renormalized one-particle retarded (advanced) Green's functions 
averaged over disorder 

(1.27) 

where T,P(E) = (T, A (E))* is the self-energy, and (ii) the random potential correlator 
(V 2 ). The self-energy E is a non-critical quantity and can be approximated by 
Hk(E) ~ —i/2r, where 1/r is the momentum relaxation rate entering the Drudc 
formula of the conductivity and isotropic scattering is assumed. 

The quantity of central interest here is the diffusion coefficient D. It follows 
from very general considerations 35 that the density-response function describing 
the change in density caused by an external space and time dependent chemical 
potential is given by 



G«- A (E) 



E-k 2 /2m 



^• A (E) 
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where D(q,ui) is a generalized diffusion coefficient. The static susceptibility (which 
is non-critical in the model of non-interacting particles) is given by xo = Np , where 
Np is the density of states at the Fermi level. The form of \ is dictated by particle 
number conservation and may be expressed in terms of G R ' A as 

x(g,w) = V$ fc fc'(g,w) + xo- (1.29) 

The two-particle quantity 

® kk '{q,u) = {G* +tk ,Gt, kL ) , (1-30) 

where G R '£ are non-averaged single-particle Green's functions, k± = (fe ± q/2, 
E±lj/2), and the angular brackets denote averaging over disorder, may be written 
in terms of the irreducible vertex function U as 

®kk>(q,w) = G R + G A _ 5fc, fc / +^J7fcfc»(q,w)$ fc "fc'(9,w) . (1.31) 

fe" 

In a diagrammatic formulation the vertex function U is given by the sum of all 
particle-hole irreducible diagrams of the four-point vertex function. By expressing 
G R G A as 

°*+ G *- = w-fe-q/m-AE,' (L32) 

where AGk = G R — G A _ and AE^ = T, R + — T, A _ one may rewrite eq. (1.31) in the 
form of a kinetic equation 



A£ fc )$ fefe , =-AG fc 



<5fcfc' + Ukk"&k"k> 



(1.33) 



fe" 

By summing cq. (1.33) over fe, fe' one finds the continuity equation 

LJ$(q,u>)-q$j(q,u)) = 2mNF (1.34) 
with the density-relaxation function 

and the current-density relaxation function 

$,(9,0;) = V^$ fefe ,(g,w), (1-36) 

TO 

fc.fc' 

where g = g/ | g |. Here the Ward identity AS fc = ^ fe , U kk 'AG k ' = U k > k AG k ' 
has been used. 33 Since the Ward identity plays a central role in the derivation of 
the self-consistent equation, we provide a short proof which does not rely on the 
perturbation expansion employed in Rcf. 33 Instead the proof follows the derivation 
of a similar Ward identity in the case of wave propagation in disordered media. 36 
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Starting from the equations of motion of the single particle Green's function before 
impurity averaging 

E + | + i0 + ^-V 2 ri - V{n)] G R (r u r 2 ; E + |) = 5{n - r 2 ), (1.37) 



E---i0 + -Vl 3 -V(r 3 ) 



G A (r 3 ,r 4 ;E--) = S(r 3 -r 4 ) 



(1.38) 



we multiply the first of these equations by G A {r 3 ,r 4 ] E — ^) an d the second by 
G R (r 1: r2\E + ^) and take the difference. We now perform the limit r% — >• r% , 
upon which the terms containing the disorder potential V(ri), i — 1,3, cancels out. 
Finally, the disorder average is taken and the result is Fourier transformed into 
momentum space, with the result 



k 



k ■ q 



^kk'(q,^) 



Gi 



(1.39) 



Comparing with eq. (1.33) it is seen that the Ward identity indeed holds. 

In the hydrodynamic limit, i.e., ujt 1, qi 1, the current density is propor- 
tional to the gradient of the density, which is expressed in Fourier space by 



$j + iqD(q,Lu)$ = 0. 



(1.40) 



In fact, multiplying cq. (1.33) by k ■ q/m and summing over k and fe', one may 
derive relation (1.40) and by comparison finds 

D /D{q,u) = l-i 1 —Y^{k-q)G^GtUkw{q^)G§G A , (fe' ■ q), (1.41) 
mn * — ' + + - 

k,k' 

where r\ = 7tNf(V 2 ) = 9 J Er is the disorder parameter, and Do = ^» 2 t is the bare 
diffusion constant. 

As the Anderson transition is approached the left-hand-side of eq. (1.41) will 
diverge for q,uj — > 0, and therefore the irreducible vertex U has to diverge, too. 
The leading divergent contribution to U is given by the set of diagrams obtained by 
using the following property of the full vertex function T (the sum of all four-point 
vertex diagrams) in the presence of time-reversal symmetry: 32 ' 3 



Ffcfc'(q, uj) — T( h _ k > 



-q)/2,(k> 



-k+q)/2(k + k', 



(1.42) 



This relation follows if one twists the particle-hole (p-h) diagrams of T such that the 
lower line has its direction reversed, i.e., the diagram becomes a particle-particle 
(p-p) diagram. Now, if time-reversal symmetry holds, one may reverse the arrow 
on the lower Green's function lines if one lets k — > —A; at the same time. This 
operation transforms p-p-diagrams back into p-h diagrams, so that an identity is 
established relating each diagram of V to its transformed diagram T T , which yields 
the above relation. 
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The leading singular diagrams of T give rise to the diffusion pole 

11 

Td= 2^N^ -iu + Dq*' (L43) 

where D is the renormalized diffusion coefficient. These diagrams are of the ladder- 
type and therefore reducible. Their transformed counterparts arc. however, 
irreducible and thus contribute to U. We may therefore approximate the singular 
part of U by 

77 sin s = I \ d 44^ 

kk ' 2ttN f t 2 -iuj + D(k + k') 2 ' K ' ' 

In low-order perturbation theory U sms is given by the "maximally crossed dia- 
grams", which when summed up give a result {7 smg ' similar to cq. (1.44), with 
D replaced by the diffusion constant Dq. When U slng '° is substituted as a vertex 
correction into the conductivity diagram, the result is exactly the weak-localization 
correction discussed in Sec. 1.1.3. The structure of the kernel Ukk' has been ana- 
lyzed from a general viewpoint in Ref. The importance of the diffusion pole for 
the Anderson localization problem was discussed in Ref. 38,39 in connection with the 
derivation of mean-held theories for disordered systems in the limit of high spatial 
dimensions. 



1.4. Self-consistent theory of Anderson localization 

It follows from eq. (1.41) that for d < 2 even the lowest-order correction in the dis- 
order parameter r/ to the inverse diffusion coefficient (obtained by replacing Ukk' by 
JJ slns '°) yields a contribution which, in principle, diverges in the limit ui — > 0. This 
infrared divergence depends crucially on the dimension d and leads to a breakdown 
of perturbation theory in dimensions d < 2. In higher dimensions the divergence 
takes place at finite disorder strength. Since the fundamental reason for the di- 
vergence of Do/D(0, 0), eq. (1.41), is the presence of diffusion poles in the kernel 
Ukk 1 , and since these diffusion poles depend on the renormalized diffusion coeffi- 
cient, Vollhardt and Wolfle 32 ' 33 interpreted eq. (1.41) as a self-consistent equation 
for the diffusion coefficient. 

By construction eq. (1.41) is in agreement with perturbation theory. An earlier 
attempt to set up a self-consistent equation in the spirit of mode-mode coupling 
theory 40 failed to reproduce the weak localization results, as it did not account for 
quantum interference effects. A later ad hoc modification of the latter theory led to 
a self-consistency scheme 41,42 which is in partial agreement with the one presented 
here, the main difference being that an additional classical (i.e., not interference 
related) mechanism of localization is included. 

When (7 sing from eq. (1.44) is substituted for U, eq. (1.41) for the diffusion 
coefficient D(u>) (i.e., in the limit q — > 0) leads to the following self-consistent 
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equation for the frequency-dependent diffusion coefficient D(lu) 



.32,33 



Do 1 , k 2 F d f 1 " _ Q d ^ 



1 + ^— dQ . 2 n( w (1-45) 

Here we assumed that a finite limit lim q ^,o D(q,u) = D(uj) exists, and that Q is 
limited to l/£ in the diffusive regime. 
Eq. (1.45) may be re-expressed as 

D{UJ) 1 - V dk 2 F - d f' 1 dQ . .ff" 1 , rxl . (1.46) 



Do ' b J ^ -ilo/D{uj) + Q1 

1.4.1. Results of the self-consistent theory of Anderson localization 

In d = 3 eq. (3.14) has a solution in the limit u — > up to a critical disorder 
strength r] c 

D = D (1--), r?<r ?c = ^=, (1.47) 

?7c V37T 



which implies the critical exponent of the conductivity s = 1. The w-dependence of 
D(u>) at the critical point is obtained as 43 

D(u) = Doiior) 1 / 3 , r) = r) c , (1.48) 

implying a dynamical critical exponent z = 3 in agreement with the exact result of 
Wegner. 19 

At stronger disorder, rj > r] c , all states are found to be localized. The localization 
length £, defined by £~ 2 = lim^^o^ i^/D^)), is found as 



1 - — 



(1.49) 



i.e., the exponent is v = 1. For general <i in the interval 2 < rf < 4 one finds Wegner 
scaling, s = i/(c£ — 2). An extension of the self-consistent theory with respect to the 
momentum dependence of the renormalizcd diffusion coefficient near the Anderson 
transition has been proposed in Ref. 44 It leads to a modified critical exponent of the 
localization exponent, v — \j(d — 2) + 1/2, which is in much better agreement with 
numerical results in d = 3. The conductivity exponent is found to be unchanged 
(s = 1), i.e., Wegner scaling is no longer obeyed. 

In dimensions d < 2, there is no metallic- type solution. The localization length 
is found as 



exp 1 

n 



1/2 



C = ci£, d = l (1.50) 

where the coefficient ci w 2.6, while the exact result is Ci = 4. 11 

The /3-function has been derived from the self-consistent equation for the length- 
dependent diffusion coefficient, where a lower cutoff l/L has been applied to the 
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Fig. 1.1. Renormalization group /3-function in dimensions d = 1, 2, 3 for the orthogonal ensemble, 
as obtained from the self-consistent theory. 45 

Q-integral in eq. (1.46). The result 45 for d = 3 dimensions in the metallic regime is 
given by 

/?(<?) = g>g c = \, (1.51) 

g tt z 

and in the localized regime by 

K9) = 1 - ^r\^z- X - TT^- . 9<9c (1-52) 
ir z g 1 + ar 1 + x 

Here x = x(g) is the inverse function of 

g = 4j(l + x)e~ :!: (l - xarctani). (1.53) 

The /3-functions in d = 1, 2, 3 obtained in this way are shown in Fig. 1.1. 

The phase boundary separating localized and extended states in a disordered 
three-dimensional system may be determined approximately by a variety of meth- 
ods. For electrons on a cubic lattice with nearest-neighbor hopping and one or- 
bital per site with random energy chosen from a box distribution in the interval 
[—W/2,W/2], the phase diagram has been determined by numerical simulations 46 
as shown in Fig. 1.2. 

Also shown is the result of an analytic expression obtained from the self- 
consistent theory 47 applied to a tight-binding model, where the Coherent Potential 
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Fig. 1.2. Phase diagram showing metallic (M) and insulating (I) regions of the tight-binding 
model with site-diagonal disorder (box distribution of width W). Dots: numerical study; 4 solid 
line: self-consistent theory. 47 The remaining lines are bounds on the energy spectrum; see Ref. 34 

Approximation (CPA) was used to evaluate the single-particle properties; no ad- 
justable parameters enter. The agreement is seen to be very good. 

1.5. Applications of the self-consistent theory of Anderson localiza- 
tion 

The self-consistent theory of Anderson localization proposed by us in 1980 32,33 was 
applied and extended to account for many of the salient features of disordered 
systems. Here we briefly review the more recent developments, not yet described in 
our review. 34 While initially the main interest had focussed on disordered electronic 
systems, in recent years the interest shifted to localization of classical waves and 
even more recently, to ultracold atom systems. We first review an extension of the 
self-consistent theory to the case of weak applied magnetic and electric fields. 



1.5.1. Effect of static magnetic and electric fields 
1.5.1.1. Magnetic fields 

One of the limitations of the self-consistent theory has been the difficulty to treat 
scale dependent contributions to the conductivity in the presence of a magnetic 
field in perturbation theory. As explained above, a magnetic field induces a "mass" 
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in the Cooperon propagator and therefore removes the localizing interference effect 
leading to localization of all states in d < 2 dimensions. On the other hand studies 
of the nonlinear a- model show that in higher (two-loop) order scale dependent terms 
appear which are generated solely by diffusion propagators (diffusons). There is, 
however, a general theorem of perturbation theory, related to gauge invariance, 
stating that the singular contribution of any diagram with one diffuson and an 
arbitrary decoration with additional impurity lines cancels within a group of related 
diagrams. 33 The way out of this apparent contradiction has only been found very 
recently. 48 

A satisfactory generalized self-consistent theory for the case of unitary symme- 
try, including the two-loop and higher contributions has not been formulated yet. 
Nonetheless there is a parameter regime of weak magnetic field B (u c = eB/mc ) 
and moderately strong disorder, lo c t << I/eft ^ 1 in which the one- loop contribu- 
tions still dominate over the two-loop contributions and a generalized self-consistent 
theory may be formulated. The most complete discussion of this approach was given 
by Bryksin and Kleinert, 49 who proposed a set of two coupled self-consistent equa- 
tions for the diffusion coefficients D p h in the particle-hole channel (diffuson) and 
D pp in the particle-particle channel (Cooperon) of a 2d system: 



i_ 1 ., w | +< ^ + _4_ ) .^ + _^_ )] , „ 54) 
^ = 1 _ 9M1+ 2^ ) . (1 . 65) 

Here i^{z) is the digamma function, g = 2/(irkp£) is the coupling constant, 
£b = (c/eB) 1 / 2 is the magnetic length, and n = l/(\/2£). The solution of these 
equations allows one to extend the results of weak localization theory, e.g., for the 
negative magnetoresistance, to the regime of moderately strong disorder, leading to 
renormalized values of the parameters of weak localization theory. Good agreement 
has been found with experimental data in that range. 49 



1.5.1.2. Electric fields 

An applied static electric field E affects the localization physics in the following 
way: electrons drifting under the influence of E experience a reduced probability of 
return, weakening the localization effect provided by interference of return paths. 
This effect is incorporated into the Cooperon dynamics, leading to a new term in 
the diffusion pole 

(1.56) 



2itN f t 2 -iu + Dq 2 + ifj, d q ■ E ' 

where fi^ = e/ (mr) is the mobility. The electric field term leads to the appearance 
of a localization transition even in dimension d = 1, 2 . Near the transition in d = 1 
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the diffusion coefficient is found as 



49 



D(E) = 



D (l-E /E), forE>E 
for E < E, 



(1.57) 



in agreement with the exact result in Ref., 50 where Eq — (2n/(eirN F Do) is the 
threshold field. In d = 2 dimensions the behavior above threshold is logarithmic: 



where E = (4ep/eK,) exp(— Trkpl/2). The relaxation of the charge current following 
a sudden switch on of the electric field has been considered in Ref. 51 There it was 



singular behavior of the Cooperon pole. 

The way in which electric and magnetic fields affect transport near the localiza- 
tion transition in anisotropic systems was studied in Refs. ' 53 

1.5.2. Anisotropic systems, films and wires 

The question of how the scaling properties of the conductance are modified in 
anisotropic systems was first addressed in Ref. 54 There it was established that 
even in the presence of an anisotropic electronic band structure and an anisotropic 
impurity-scattering cross section the one-parameter scaling theory holds. The ratios 
of the components of the conductivity tensor are invariant under scaling, implying 
that the geometric mean of the conductivity components plays the role of the scaling 
quantity. This feature is preserved by the self-consistent theory. Numerical studies 
of anisotropic systems 55,56 appeared to cast doubts on the one-parameter scaling 
hypothesis. However, a later more careful study of the problem in d = 2 dimensions 
showed that indeed one-parameter scaling is obeyed: 57 the ratio of the localization 
lengths (in the direction of the principal axes) turns out to be proportional to the 
square root of the ratio of the conductivities. A comparison with the self-consistent 
theory in the somewhat simpler form of the "potential- well analogy" 58 showed again 
qualitative agreement. Localization in anisotropic systems has also been considered 
in a model with anisotropic random potential correlations, and the phase diagram 
has been mapped out within an extension of the self-consistent theory. 59 The 
same authors explored the consequences of finite-range correlations of the random 
potential within a generalization of the self-consistent theory. 60 

A somewhat different but related question is the behavior of the conductance 
of a film of finite thickness, or a wire of finite diameter. There is no doubt that 
in the thermodynamic limit these systems behave like true 2d or Id systems. It 
is, however, interesting to understand how this behavior is approached. Numerical 
studies of metallic disordered films as a function of film thickness seemed to indicate 
a localization transition as a function of thickness, 61 in contradiction to the results 
of the self-consistent theory applied to this system. A further study by the same 




(1.58) 



found that the current has a long time tail cc t x l 2 as a consequence of the infrared 
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authors 62 on systems of finite thickness in a magnetic field explored the possibility of 
a derealization transition controlled by both the thickness and the magnetic field. 
The transitions obtained are pseudo-transitions marking a crossover from strong to 
weak localization, as confirmed in a later more accurate numerical study. 63 

1.5.3. Anderson localization of classical waves 

The concept of the self-consistent theory of localization can be carried over to 
the case of propagation of classical waves in disordered media. Here we sketch 
the formulation following the presentation of Kroha, Soukoulis, and Wolflc. 64 For 
scalar waves propagating in a medium of randomly positioned point scatterers of 
density nj — a -3 , modelled by spheres of volume V s , the average phase velocity 
is given by c p h = cq[1 + (V s /a 3 )Ae]~ 1 / 2 , where cq is the bare phase velocity and 
Ae characterizes the strength of the scattering ( "dielectric contrast" ) . The Green's 
function of the wave equation is defined as 

G fe (w) = [Go 1 M-S fe (cj)] ! (1.59) 

where Gq 1 (uj) = uj 2 — c 2 k 2 . The self-energy £ may be determined within the Coher- 
ent Potential Approximation (CPA) (see, e.g., Ref. 65 ) provided it is independent of 
k. Then the bare diffusion constant is found as 



Do = 2cH^G- 1 M £(ft • q) 2 (ImG£) 2 . (1.60) 



The renormalized diffusion coefficient may be shown to satisfy the self-consistency 
equation 

nm\ n or I x c o l2 ImS D (n) V""7i ~\ ImG fc (ImG fc >) 2 

m = Do (ImGo)2 Do E( fc -g) _^ +jW+fc , )2 ( fc ■«)• 

(1.61) 

Here fl is the external frequency while uj is the frequency of the waves which enter 
in one-particle quantities. This equation can be solved in the limit Q — > to 
obtain the diffusion coefficient in the delocalized phase and the localization length 
£ = \imQ^a[D(fl) / (—ifl)] 1 ^ 2 in the localized phase. One finds that it is much harder 
to localize classical waves as compared to electrons, and there is only a narrow region 
of the phase diagram (at reasonable contrast Ae ) where localization is found. 64 

A more realistic theory of the propagation of electromagnetic waves in disordered 
materials with loss or gain mechanisms keeping the vector character of the fields 
has been worked out by Lubatsch, Kroha, and Busch. 66 We briefly sketch the 
main results here. The electric field amplitude E^{r) of an electromagnetic wave of 
frequency a; in a medium with random dielectric constant e(r; uj) = e(u>) + Ae(r; uj), 
obeys the wave equation 

Vx(VxE u (r))-^e(r;uj)E u (r) = ujJ u (r), (1.62) 
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where (Ae(r;w)) = 0. In the following the random part of the dielectric function 
will be modelled as Ae(r;w) = — (c 2 /uj 2 )h(uj)V(r). The Green's functions of the 
wave equation, after disorder averaging, are defined as 

,:2 



Gf\u) = [{^e-k")P-^ A ^)]-\ (1.63) 



Here G and the self-energy X are (3 x 3) tensors and P = 1 — k<Eik is the projector 
onto the transverse subspace (here and in the following the hat symbol denotes 
a unit vector). Transport properties are contained in the two-particle correlation 
function (a tensor of rank four) 

**fc/(g,fi) = (G« x+ ® Gt tk >_) > (i- 64 ) 
where k± = (fe ± q/2,u) ± f2/2) , etc., which obeys the Bethe-Salpeter equation 

* fcfe , (g, fi) = Gf + ® [<W + J2 u *k" («> ")*fc»fc' («= 0)] ■ (1-65) 

fe" 

As in the case of electrons in a random potential considered above, the Bcthc- 
Salpctcr equation may be converted into a kinetic equation for the integrated in- 
tensity correlation tensor 4>fc(q, fi) = J2k' ^kk'{q, of the form 

(AG^(w)-AE fc )* fc = - AG fc [l® l + ^C7 fcfc »* fc »], (1.66) 



where AG^ (u) = [G^]" 1 ® 1 - 1 ® [G^]" 1 , AS fc = E« $ 1 - 1 $ £ 



and AG fc = Gf ® 1 - 1 ® G^_ . 



The kinetic equation serves to derive the energy conservation equation and the 
equivalent of Fick's law: 

Ti("J 

J E (q,Q)=iP E (q,Q)D(n)-q. (1.68) 

Here 

P £ ( 9 ,r!) = ( W / Cp ) 2 ^* fc ( q ,r!) (1.69) 
fe 

is the energy-density relaxation function, with c p as the renormalized phase velocity, 
and 

£(fc.g)* fc (g,n) (1.70) 
fe 

is the energy-current density relaxation function, with v E {uS) as the energy trans- 
port velocity; for the definitions of c p and v e (oj) we refer the reader to Ref. 66 
When energy absorption by the medium is taken into account (as expressed by the 
imaginary part of the dielectric function), or conversely, if a medium with gain is 
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considered, energy is not conserved, as expressed by the loss/gain rate ■ The 
energy-diffusion coefficient tensor D(Cl) is found as 

D(Sl) = ^v e (u)It, (1.71) 
where the tensor of transport mean free path is given by 

It = ^(A + k)1. (1.72) 
Here the main contribution to I has a form which is analogous to cq. (1.41): 

r 1 =a- 1 J2(k-q)AG k U kk >>(0 i n)AG k >(k' -q), (1.73) 

and ai and A arc defined in Ref. 66 The quantity k describes scattering caused by 
a mismatch of absorption/ gain between the scattering objects and the medium. 

The energy density propagator Pe(<7, f2) in the limit of small q, O follows from 
eqs. (1.67), (1.68) as 

p E ( g , n) = [si + + *'« ■ D ity ■ g] _1 <%, a). (1.74) 

Replacing U kk " (0,Q) by its singular part proportional to the diffusion propagator 
Pe(q, fi) one arrives at a self-consistent equation for the diffusion coefficient ten- 
sor. The latter provides a framework for the description of the interplay between 
localization and stimulated emission in materials with gain, i.e., the problem of the 
random laser. 67 

The predictions of the self-consistent theory have also been probed by compar- 
ison with numerical results 68 for transmission of waves in unbounded Id and 2d 
systems and through strips of finite width. Good overall agreement is found. 

The localization of phonons and the ultrasound attenuation in layered crystals 
with random impurities has been studied within the self-consistent theory in Ref. 69 

1.5.4. Transport through open interfaces 

Most of the discussion so far considered transport in infinitely extended systems, 
with the exception of the scaling theory for systems of length L. In some cases, how- 
ever, transport through plate-shaped systems in the direction perpendicular to the 
plate surface is of interest. As pointed out by van Tiggelen and collaborators, 70,71 
the weak localization physics changes near an open boundary, as the finite prob- 
ability of escape through the interface diminishes the return probability necessary 
for interference. In the framework of the self-consistent theory this effect may be 
taken into account quantitatively. To this end it is useful to express eq. (1-45) in 
position-energy space as 

l + 27r^— C(r,r'), (1.75) 



D(lu) m 
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where C(r,r') is a solution of the diffusion equation, and a cut-off Q < l/£ to the 
momentum was applied to the spectrum of the Q-modcs in eq. (1.45) 



The above formulation now allows one to describe position-dependent diffusion pro- 
cesses, as they appear near the sample surface in a confined geometry, e.g., trans- 
mission through a slab. In that case the diffusion coefficient may be assumed to 
be position dependent, D = D(r,ui). Then C(r,r/) obeys the modified diffusion 
equation 72 



The solution is subject to an appropriate boundary condition at the surface of the 
sample. A microscopic derivation of the above equation in diagrammatic language 
was given in Ref. 73 Further confirmation of the theory was obtained in, 74 where 
the above equations was derived within the nonlinear er-model framework in the 
weak coupling limit. The theory accounts very well for the localization properties 
of accoustic waves transmitted through a strongly scattering plate. 72 

It is natural to ask whether a position dependent diffusion coefficient will change 
the critical behavior obtained from the scaling properties of the conductance of 
finite size samples. This question was addressed in Ref. 75 with the result that 
the critical exponents are unchanged and the /3-function is hardly modified by the 
improved approximation. The scaling of the transmission coefficient for classical 
waves through a disordered madium near the Anderson transition was considered 
within the position dependent self-consistent theory in Ref. 76 

The transmission of microwave pulses through quasi-onc-dimcnsional samples 
has been measured recently and was analyzed in terms of the self-consistent the- 
ory. 77 It was found that while the self-consistent theory can account very well for 
the propagation at intermediate times, it fails at longer times when the transport 
occurs by hopping between localized regions. 

Anderson localization of atoms in a Bosc-Einstcin condensate released from a 
trap and subject to a random potential has been considered in the framework of the 
self-consistent theory in Ref. 78 The authors show that the scaling properties govern 
the dynamical behavior of the expanding atom cloud, so that the critical exponents 
determine the power law in time obeyed by the expanding cloud size. 

1.6. Conclusion 

Anderson localization in disordered systems continues to be a very lively field of 
research. Current investigations do not concentrate so much on disordered elec- 
trons but on classical waves (light, electromagnetic microwaves, acoustic waves), or 
ultracold atoms in the presence of disorder. Although the fundamental concepts of 
Anderson localization are well understood by now, there still remain a number of 



\-iu) + L»(w)V 2 ]C(r, r') = S(r - r'). 



(1.76) 



+ VD(r,u)V]C(r,r') = 5(r-r'). 



(1.77) 
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open questions. Some of them are related to the analytical theory of critical prop- 
erties near the Anderson transition. Others concern the quantitative description of 
realistic materials, e.g., the question under which conditions light or acoustic waves 
become localized. The self-consistent theory of Anderson localization has been, and 
will continue to be, a versatile tool for the investigation of these problems. It allows 
one to incorporate the detailed characteristics of the system such as the energy 
dispersion relation, the particular form of disorder, the shape of the sample, and 
loss or gain mechanisms in an efficient way. The self-consistent theory is not only 
applicable to stationary transport problems, but also to dynamical situations such 
as pulse propagation or the behavior after a sudden switch-on. 

As Anderson localization is a wave-interference phenomenon, the limitations 
of phase coherence are an important subject of study in this context. By now 
Anderson localization has been observed in many different systems beyond doubt. 
On the other hand, the observation of the Anderson transition itself is a much more 
challenging task. Here the recent investigations of classical waves and atomic matter 
waves offer fascinating, new perspectives which will undoubtedly lead to a deeper 
understanding of the localization phenomenon. 
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